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ABSTRACT 

We investigate a population of transiting planets that receive relatively mod- 
est stellar insolation, indicating equilibrium temperatures < 1000 K, and for 
which the heating mechanism that inflates hot Jupiters does not appear to be 
significantly active. We use structural evolution models to infer the amount of 
heavy elements within each of these planets. There is a correlation between the 
stellar metallicity and the mass of heavy elements in its transiting planet (s). 
It appears that all giant planets posses a minimum of ~ 10-15 Earth masses of 
heavy elements, with planets around metal-rich stars having larger heavy element 
masses. There is also an inverse relationship between the mass of the planet and 
the metal enrichment [Zpi/Zstar), which appears to have little dependency on the 
metallicity of the star. Saturn- and Jupiter-like enrichments above solar compo- 
sition are a hallmark of all the gas giants in the sample, even planets of several 
Jupiter masses. These relationships provides an important constraint on planet 
formation, and suggests large amounts of heavy elements within planetary H/He 
envelopes. We suggest that the observed correlation can soon also be applied to 
inflated planets, such that the interior heavy element abundance of these planets 
could be estimated, yielding better constraints no their interior energy sources. 
We point to future directions for planetary population synthesis models and sug- 
gest future correlations. This appears to be the first evidence that extrasolar 
giant planets, as a class, are enhanced in heavy elements. 

Subject headings: planetary systems 



1. Introduction 

Transiting exoplanets are valuable for planetary characterization because they allow us 
to measure their masses through stellar radial velocity or other dynamical measurements. 
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as well as their radii from the transit light curve. Together, these yield a planet's bulk 
density. In principle, this information could be used to determine a planet's composi- 
tion as increasing the mass fraction of heavy elements increases the density. This appar- 
ently straightforward method has been difficult to implement, however. Transit observa- 
tions have revealed that most of the highly irradiated "hot Jup iters" are inflated to large 
radii beyond what is expected from simple models. The reason for this effect has not been 



determined; a variety of adc 
anisms have been proposed ( 


itional internal ener 


gy sources or contraction-stalling mech- 


Guillot & Showman 


2002 




Jackson et al. 


2009: 


Batvein et al. 


2011; 


Chabrier & Baraffe 


2OO7I: 


Arras & SocratesI 


2010f 


. 


Since an inflated radius decreases 



a planet's density, the heating mechanism acts in opposition to the effect of adding heavy 
elements to the planet. Therefore interior composition for a transiting planet i s generally 
left e r itirely unknown, un l ess planets are found to be dramatically overdense (e.g. 
2OO5I : iFortnev et aPbood : lleconte et aPboOoh . 



Sato et al. 



For the over-inflated planets it is possible to flnd a relation between the heavy elements 
in the planets and the metallicity of the stars by making an assumption about the relation- 
ship between the in cident stellar radiation and the unknown power input into the planet 
( iGuillot et al.l 120061 ). A similar relationship has b een found by using ad hoc enhanced at- 
mospheric opacity to slow planetary contraction ( iBurrows et al.l 120071 ). These studies are 
intriguing, although the resulting planet-star metallicity relationship is dependent on the as- 
sumed behavior of the unknown radius inflation mechanism. Since it is well-known that our 
solar system's four giant planets possess at least 10 — 15 of heav y elements within their 
interi ors, making them enriched compared to the Sun's composition (IFortney &: Nettelmann 
2OIOI ). it is paramount to determine the composition of giant exoplanets to understand the 
structure and formation of these planets as a class of astrophysical objects. 

Empiri cally, the unknown h eating mechan ism affects the close-i n planets at high incident 
stellar flux ( iKovacs et al.ll2010al ) or planet Teg ( iLaughlin et al.ll201ll ) . This is shown in Figure 
1, where we plot the observed planet radii as a function of their average incident stellar flux. 
A thermal evolution model for a 1 Jupiter- mass (1 Mj) planet with no extra heating source 
other than the effect of incident radiation on the planet's atmosphere is plotted with a h eavy 



element core (25 Earth masses, M^, dotted) and without (solid) ( iMiller et al. 



the sample is divided into two regions l|: those 



20091). At 



approximately {F) = 2 x 10 erg s~ cm"', tne sample is divided into two regions LJ: tnose m 
the higher flux region, most of which require extra heating to explain their radii, and those in 
the lower flux region, in which no inflated planets are found. This empirically suggests that 



^{F) = 2 X 10^ erg s"^ cni"^ corresponds to an equilibrium temperature of 990 K for an Bond albedo of 
0.1 and efficient heat transport between the day and night side. This temperat ure is quite similar to that 
for which Ohmic heating is suggested to become important (jBatygin et al.ll201ll ) 
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Fig. 1. — Planet radius as a function of average incident stellar flux. Planets are colored 
according to their mass. Model planet radii are plotte d for a 1 Mj plane t at 4.5 Gyr 



witho ut a core (solid) and with a 25 core (dotted) (IFortney et al. 



2007: Miller et al. 



20091 ). Although the extra heating source is not well-determined, it is clear that it is more 



important at larger incident fluxes. We choose a cutoff of (F) < 2 x 10* erg s ^ cm ^ in 
order to obtain the largest sample of non-infiated planets. This corresponds to a planetary 

Teq < 1000 K. 



the heating mechanism does not significantly contribute to the energy budget at low incident 
flux. Therefore, for this sample of 14 transiting exoplanets, we can neglect the heating 
mechanism and use our structural evolution models to estimate a planet's composition. 




Fig. 2. — The stellar metallicity and inferred planet heavy element mass for exoplanets 
within our incident flux cut. The required heavy elements are from the "Average Case" in 
Table 1. (See text.) Planets are numbered corresponding to the entries in Table 1. Th e 
rarity of gas giants around metal-poor stars is well established (IFischer fc Valentil l2005l ). 
Using a least-squares fit, we find the relation logM^ = (0.82 ± 0.08) + (3.40 ± 0.39) [Fe/H] 
and a reduced Chi-squared value of 1.95. The fit excludes HAT-P-12b (planet 13) and 
includes Jupiter and Saturn. However, we do not expect this relation to hold at the lowest 
metallicities, where it may become fiat at ~ 10-15 M,^. 



2. Model and Method 

We consider two limiting types of planetary structures. We consider planets where all 
of the heavy elements are in an inert core with an adiabatic solar metallicity convective 
envelope above (layered model). We also consider a structure where heavy elements are 
uniformly mixed with the hydrogen and helium and the planet is fully convective (mixed 
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model). The primary heavy elei nent composition is a mixture of 50% rock and 50% ice using 
the equation of state ANEOS (IThompsonl Il990l ) . By considering the two extreme cases of 
having all of the heavy element masses in the core or envelope, we bracket possible interior 
models of giant planets. For Jupiter, models that match gravity field constraints generally 
find t hat most of its heavy eleme nts are in the envelope while for Saturn most are in the 
core ( iFortney fc Nettelmannll2010l ). 



A com plete description o f the thermal evolution model can be found in iFortney et al. 



(2007 


) and 


Miller et al. 


(2009) 



an inert core, 2) an adiabatic convective envelope (wh ere heavy elernents m ay be mixed in) 



and 3) a solar-metallicity non-grey atmosphere model (IFortney et al.ll2007( ) that includes the 
atmospheric extension to the transit radius. The primary effect of heavy elements either in 
the core or in the convective envelope is mainly to decrease the planet's radius at every time. 

For each planet, the amount of heavy elements is determined under the constraint that 
the predicted model transit radius agrees with the observed radius at the observed age and 
incident flux. The average incident flux that a planet receives is given by 

~ w7r^ 

where is the luminosity of the star, a is the semi-major axis of the orbit and e is the 
eccentricity of the orbit. This analysis was performed on all planets that met our average 
incident flux cut {F) < 2 x 10* erg s~^ cm~^ and had a mass greater than 20 M^-since our 
model is primarily designed to describe giants with masses greater than Neptune. 

Note these heavy element masses should be taken as minimum masses since if the planet 
is internally heate d or if higher atmospheric opacities (du e to metal-enhanced atmospheres) 
slow the cooling (llkoma et al.l l2006l : iBurrows et al.l 120071 ). then a planet would have more 
heavy elements than found here. 

The required heavy element mass to fit the radius is determined as the average of the 
layered and mixed cases. Each of the the observed system parameters (i?p, age, a, Mp) has 
an associated error on its published value. The propagated error on the heavy element mass 
{an) is given by: 

2 



H 





2 




2 




dRp 




dAge 


O^Agc + 


da 



+ 



a 



Mr 



Mr - Mp. 



(2) 



where an^, crAgc, <^a-, and om^ are the observationally determined errors in planet radius, 
system age, semi-major axis, and planet mass respectively. The derivatives (calculated 
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at the observed planet parameters assuming core heavy elements) describe the sensitivity of 
the predicted heavy element mass with respect to changes in a given parameter, X. The 
final term of the expression is the uncertainty due to the unknown structure of the planet. 
Mc and Mgnv are the predicted heavy element masses if the heavy elements are within the 
core, or the envelope, respectively. 

We use the metallicity of the star [Fe/H] as given in each paper in Table 1. For each 
system, we compute the heavy element mass fraction Z^tax = 0.0142 x 10^^*^/^^ - assuming 
that the total heavy element compositio n of other systems sc ales with their iron abundance, 
normalized to the solar metalicity as in lAsplund et al.l (120091 ). 



3. Findings 



In Figure 2 we plot the stellar metallicity, [Fe/H], against the planet heavy element mass 
for each of these systems. Using a least squares fit, we find that logM^ = (0.82 ± 0.08) + 
(3.40 ± 0.39) [Fe/H] for stars with [Fe/H] > -0.05. The reduced Chi-squared value of 1.95 
implies that not all of the scatter can be explained by observational error. We expect a fairly 
flat relation (the dotted line in Figure 2) at subsolar stellar metallicity if 10-15 of heavy 
elements are needed to trigger planet formation. In Table 1 we list the planets and observed 
parameters used. For each planet, we list the average predicted heavy elements between the 
core model and mixed model with the 50-50 rock-ice composition. 

We have examined the sensitivity of our flndings to alternate choices for the heavy 
element EO S, and the differe nces are small. On the low density E OS end, we have used 
100% wat er (Thompson 1 199 ol), a nd on the high density end 2/3 rock ( Thompson 1990) and 
1/3 iron (ILyon &: Johnsonlll992l ). These generally lead to Zpi/Zgtar that are 10-20% larger, 
for the water EOS, and 10-20% smaller, for the rock/iron EOS, than those found in the last 
column of Table 1. For example, for the HAT-P-17b average case, the best-fit heavy element 
mass of 16.0 ± 7.7 increases to 17.7 ± 4.8 for pure water, and 14.1 ± 3.9 for 
rock/iron. This yields Zpi/Zgtar values of 7.4 ± 2.5 and 5.9 ± 2.0 for the water and rock/iron 
cases, respectively, very similar to the 6.7 ± 3.5 value in Table 1. 

Perhaps the clearest way of looking at this sample is to compare the planet mass against 
the inferred heavy element mass or heavy element enrichment {Z^i/ Zstar) , as shown in Figure 
3. The mass of heavy elements appears to increase with planet mass. On the other hand, 
the heavy element enrichment decreases with increas ing planet mass consistent w ith the 
pattern found in the solar system's four giant planets (IFortney fc NettelmannI 120101 ). All of 
these planets are consistent with being enriched in heavy elements and many of these must 
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be enriched significantly. 

Within the cluster of Saturn-like planets around 0.2 Mj, the most heavy element-rich 
planet (10, CoRoT-8b) orbits the most metal-rich parent star. For the planets more massive 
than Jupiter, the planets harbor large amounts of heavy elements, and orbit around metal- 
rich stars, which explains some of the steep slope from Figure 2. It is interesting to note 
that the relationship of Zpi/Zgtar as a function of planet mass, appears fairly independent of 
stellar [Fe/H]. 
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Fig. 3. — Top: The planet mass and heavy element mass for our sample. Planets are 
colored by metallicity in three bins: [Fe/H] < 0.0 (red), < [ Fe/H] < 0.2 (green), and 
[Fe/H] > 0.2 (blue). Jupiter and Saturn are also shown in black (jGuillotlll999l ). This plot is 
consistent with a minimum heavy element mass of 10-15 M^, with increasing heavy element 
masses for larger mass planets. Bottom: The planet mass and heavy element enrichment 
ratio Zpi/Zstar- Lower mass planets are more metal enric hed, but have less total hea vy 
elements, which is consistent with the solar system's giants (IFortney fc Nettelmannll2010l ). 



Table 1: Table of Planets with low incident flux 



Number 


Name 




Mass 






Rad 


us 


Age 




< 


F) 


[Fe/H] 


Core Model 


Mixed Model 


Average Model 








-^pl/ ^star 


References 








(Mj) 






(«j) 




Gyr 


erg 


s 


1 cm-2 




(Me) 


(A-%) 


(Me) 


















1 


HD 80606 b 


3 


940± 





110 


1 


030± 





036 


7.o±li! 


1 


67 


XlO' 


0.390± 0.030 


87.0± 62.6 


84.5± 62.6 


85. 8± 62.6 





068± 





050 


2.0± 1.4 


[A], [A'] 


2 


CoRoT-9 b 





840± 





070 


1 


OSOi 





040 


4.0±y 


6 


58 


XlO'' 


-O.OlOi 0.060 


11. 1± 17.6 


10. 3± 17.6 


10. 7± 17.6 





040± 





066 


2.9± 4.8 


[B] 


3 


HD 17156 b 


3 


212± 





007 


1 


087± 





006 


3.4±8:| 

3.0±l;° 


1 


96 


XlO* 


0.140± 0.080 


38. 4± 9.3 


37. 8± 9.3 


38. 1± 9.3 





037± 





009 


1.9± 0.6 


[C] 


4 


Kepler-9 b 





252± 





013 





842± 





069 


8 


11 


XlO'^ 


0.120± 0.040 


31. 0± 9.4 


25. 0± 9.4 


28. 0± 9.9 





349± 





124 


18. 7± 6.9 


[D] 


5 


Kepler-9 c 





171± 





013 





823± 





067 


3 


14 


XlO^ 


0.120± 0.040 


20. 6± 6.4 


16. 5± 6.4 


18. 5± 6.7 





341± 





127 


18. 2± 7.0 


[D] 


6 


CoRoT-10 b 


2 


750± 





160 





970± 





070 


2.o±J:« 


5 


38 


XlO'' 


0.260± 0.070 


192. 0± 93.8 


172. 3± 93.8 


182. 1± 94.4 





208± 





109 


8.1± 4.4 


[E] 


7 


HAT- P- 15 b 


1 


946 ± 





066 


1 


072± 





043 


6.8±?:| 


1 


51 


XlO** 


0.220± 0.080 


22. 6± 37.1 


21. 5± 37.1 


22. 0± 37.1 





036± 





060 


1.5± 2.6 


[F] 


8 


HAT- P- 17 b 





530± 





018 


1 


010± 





029 


8 


91 


XlO^ 


O.OOOi 0.080 


16. 9± 7.7 


15. 1± 7.7 


16. 0± 7.7 





095± 





046 


6.7± 3.5 


[G] 


9 


WASP-8 b 


2 


240± 





080 


1 


038± 





047 


4.o±;| 


1 


79 


XlO* 


0.170± 0.070 


76. 6± 52.8 


70. 1± 52.8 


73. 4± 52.8 





103± 





074 


4.9± 3.6 


[H]00 


10 


CoRoT-8 b 





220± 





030 





570± 





020 




1 


22 


XlO* 


0.300± 0.100 


55. 2± 8.1 


45. 2± 8.1 


50. 2± 9.5 





717± 





167 


25. 3± 8.8 


[I] 


11 


HAT- P- 18 b 





197± 





013 





995± 





052 


12A±ii 


1 


18 


XlO* 


O.lOOi 0.080 


6.8± 4.7 


5.8± 4.7 


6.3± 4.7 





100± 





076 


5.6± 4.4 


[J] 


12 


HAT- P- 11 b 





081± 





009 





422± 





014 


6.5±1:? 


1 


31 


XlO* 


0.310± 0.050 


23. 5± 2.7 


20. 4± 2.7 


22. 0± 3.1 





854± 





154 


29. 5± 6.4 


[K] 


13 


HAT- P- 12 b 





211± 





012 





959± 





030 


2 1 2.0 


1 


90 


XlO* 


-0.290± 0.050 


17.7± 4.2 


14. 5± 4.2 


16. 1± 4.5 





239± 





069 


32. 9± 10.3 


[L] 


14 


GJ 436 b 





074± 





005 





377± 





009 


6.o±l:fi 


4 


03 


XlO^ 


-0.030± 0.200 


22. 1± 1.6 


19. 5± 1.6 


20. 8± 2.1 





8S8± 





108 


67.0± 40.0 


[M] 



[1] References: [A] = iHidas et al.ll2010|). f A'1 = ]Nordstrom et al.ll2008|). fB l^ jPeeg et al.ll201()). fCl = jNutzman et al.ll201Cl). f D1 = iHolman et al.ll2010l). [E1 = iBonomo et al.ll201d). 
FF1 = iKovaes et al.l[2'010bD . [G1 = <Howard et al.ll20ia ). [Hl^ iOueloz et al.ll201CD . fI1 = iBord c et al. 201^ f J1 = iHartman et al .11 20101 ). fKl^ jBakos et al.ll20icl) . fL1 = iHartman et al.ll2009l ). [M] 

= jTorres et al.ii2008l ) 
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Discussion 



The lower-irradiation planets are important for a better understanding of giant planet 
structure and formation since they allow us to probe the composition independent of ma- 
jor assumptions. By studying the relationship between stellar metallicity, planet mass, and 
heavy element mass within a planet, we will be able to test predictions of any planet for- 
mation theory, a n d specifically aga i nst p r edictions from populati on synthesis models such as 
Ida fc LinI ( 120101 ) ; iMordasini et al.l (120091 ) ; iThommes et al.l (120081 ) . These models are now be- 
ing extended beyond just planet mass vs. period, to ex plore the relationship between [Fe/H] 
and planet heavy element mass (IMordasini et al.ll2009l ). 



The core accretion formation mechanism requires that heavy elements form a core 
of ar ound ~ 10 Mm, wh ich is followed by the accretion of gas from the protoplanetary 
disk (jPoUack et al.l 119961 ). Some models fo r Jupiter indicate that most of its heavy el- 
ements are in its envelope, not in its core (IFortney &: Nettelmanru |2010| ). Alternatively, 
gravitational instability could result in enhanced metallicity through planetesimal accretion 



(Guillot fc Gladman 



Boley fc Durisen 



2000 



2010; 



Helled et al. 



2006; Helled fc Schubert 2009; Helled et al. 2010 



NayakshinI l2010l ). The planets in our population, even the more 
massive planets, are consistent with being enhanced in heavy elements relative to their 
parent star. This enrichment in heavy elements is a distinguishing c haracteristic between 
planets and low ma ss brown dwarfs with more solar-like abundances (jChabrier et al.l 120071 : 
Leconte et all 120091 ). 



If this emerging relationship between stellar metallicity and planetary heavy elements 
continues to hold with additional data, then the relationship could be used to determine the 
amount of heavy elements in a given inflated hot Jupiter with some confidence, based only on 
the parent star metallicity and planet mass. This would be powerful as it would allow for a 
straightforward determination of the additional energy needed to explain a planet's inflated 
radius. The additional energy source could then be derived for each inflated planet, as a 
function of planet mass and irradiation level, a nd could be compared to model predictions 
JCuillot fc ShowmaiJ[2002l ; iBatvrin et aDboilh . 



As additional data accumulate, modifications to the relations presented here could be 
in order. Perhaps a spread in Zp]/Z^ tar could be due to orbita l period, which could tie into 
the planet's dynamical environment (iGuillot fc Gladmanll2000l ). These relationships may be 
interesting to analyze in systems with multiple transiting planets. Another aspect related to 
orbital evolution is possible: Perhaps differences in heavy elements could be seen between 
planets that are well-aligned or mis-aligned with their stellar spin axis, as measured by the 
Rossiter-McLaughlin Effect (iGaudi fc Winnll2007l ). as these planets may have taken different 
paths to their current orbits. The relationship to stellar mass could also be investigated. 



- 10 - 



The planet formation process and composition may also be a function of th e types 



of heavy elements that are in a protoplanetary disk. Previously, iRobinson et al.l (|2006[ ) 
showed empirically that [Si/Fe] or [Ni/Fe] are correlated with the existence of planets for a 
fixed [Fe/H]. Theoretically, they also showed that ice-rich disks tend to form cores faster. 
Therefore, as the sample of cooler planets in this domain increases, it will also be interesting 
to test how these planet composition trends are a function of [oi/H] or on [Si/H], [0/H] or 
[C/H]. It may be possible to constrain the composition of the planetary heavy elements from 
such studies. 

In closing, we find evidence from a sample of 14 transiting giant planets that these 
planets, as a class, are enhanced in heavy elements. The large heavy element abundances 
found indicate that all heavy elements cannot be found solely in a core. If the solar system 
and planet formation models are a guide, then, in addition to their dense cores, the H/He 
envelopes of these planets will be enhanced in heavy elements as well, which can be tested by 
obs ervations of the atm ospheres of planets via transit or direct imaging spectroscopy (see. 



e.g. iMarlev et al.l 120071 ). 



The trends identified here, that independent of stellar metallicity, all giant planets have 
a heavy element mass of 10 or larger, that the abundance of heavy elements in giant 
planets increases steeply with stellar metallicity, that Jupiter-like enhancement over solar 
abundances are standard for gas giants, and that more massive planets tend to have lower 
enrichment, could be enhanced or refuted by additional detections of transiting planets with 
equilibrium temperatures less than 1000 K. These longer-period systems will continue to 
be detected from the ground and recently NASA's Kepler spacecraft identified dozens of 
candidates fo r a potentially dram atically larger sample of these less- irradiated transiting 



giant planets (IBorucki et al.l 1201 ll ). 



We thank Mark Marley, Kevin Schlaufman, James Guillochon, Philip Nutzman, and 
Eliza Kempton for providing feedback and encouragement. JJF acknowledges the support 
of NSF grant AST-1010017 and an Alfred P. Sloan Research Fellowship. 
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